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Abstract Like the solar corona, the external magnetic field of magnetars is twisted 
by surface motions of the star. The twist energy is dissipated over time. We discuss 
the theory of this activity and its observational status. (1) Theory predicts that the 
magnetosphere tends to untwist in a peculiar way: a bundle of electric currents (the 
"j-bundle") is formed with a sharp boundary, which shrinks toward the magnetic 
dipole axis. Recent observations of shrinking hot spots on magnetars are consistent 
with this behavior. (2) Continual discharge fills the j-bundle with e ± plasma, main- 
taining a nonthermal corona around the neutron star. The corona outside a few stel- 
lar radii strongly interacts with the stellar radiation and forms a "radiatively locked" 
outflow with a high e ± multiplicity. The locked plasma annihilates near the apexes 
of the closed magnetic field lines. (3) New radiative-transfer simulations suggest a 
simple mechanism that shapes the observed X-ray spectrum from 0. 1 keV to 1 MeV: 
part of the thermal X-rays emitted by the neutron star are reflected from the outer 
corona and then upscattered by the inner relativistic outflow in the j-bundle, produc- 
ing a beam of hard X-rays. 



1 Introduction 

Term "magnetars" was coined for neutron stars with ultrastrong magnetic fields B ~ 
10 14 - 10 15 G (Duncan & Thompson 1992; Paczyhski 1992). There is substantial 
evidence for ultrastrong fields in soft gamma repeaters (SGRs) and anomalous X-ray 
pulsars (AXPs) (see e.g. reviews by Woods & Thompson 2006; Mereghetti 2008). 
Alternative models for SGRs and AXPs are not discussed here. 

The current magnetar catalogue contains 18 objects (13 confirmed and 5 candi- 
dates)0Host supernova remnants were identified for some of them, confirming their 
young ages t ~ 10 3 — 10 s yr. The young age, combined with the currently observed 
number of objects, implies that <; 20% of all neutron stars are born as magnetars. 
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The rotation rates of all observed magnetars are moderate by neutron-star stan- 
dards. Their spin periods P are in a rather narrow range from 2 to 12 s. On the 
other hand, their spindown rates P show large temporal variations, usually related 
to the X-ray outbursts. This activity contrasts with the traditional view of neutron 
stars as passive stellar remnants. Unlike ordinary pulsars, whose activity is associ- 
ated with the open field lines that extend to the light cylinder, the magnetar activity 
must be generated in the closed magnetosphere. Its energy output greatly exceeds 
the rotational energy of the star and must be fed by magnetic energy. 

A distinct feature of these objects is their huge nonthermal luminosity, exceeding 
the spindown power by 1-3 orders of magnitude. Between the outbursts, magnetars 
display persistent or slowly decaying X-ray emission with luminosities L ~ 10 34 — 
10 35 erg s _1 . Two peaks are observed in the X-ray spectrum: one near keV and the 
other above 100 keV (see e.g. Kuiper et al. 2008). The energy fluxes in the two 
peaks are comparable. The keV peak is dominated by emission from the hot surface 
of the neutron star. The 100-keV component is clearly nonthermal and demonstrates 
the existence of a hot corona around magnetars, which provides a large fraction of 
the total observed luminosity. For comparison, the corona of the sun radiates only 
~ 10~ 6 of the solar luminosity. 

The observed magnetospheric activity must be generated by motions of the star 
surface. The crust of a magnetar is stressed by the ultrastrong magnetic field and can 
yield to the stresses through episodic starquakes or a slow plastic flow (Thompson & 
Duncan 1996). The magnetosphere is anchored in the crust and inevitably twisted 
by the crustal motions, resembling the behavior of the solar corona. It becomes 
non-potential, V x B / 0, and threaded by electric currents (Thompson et al. 2002, 
hereafter TLK02; Beloborodov & Thompson 2007, hereafter BT07). The currents 
can flow only along B, and the twisted magnetosphere remains nearly force-free, 
j x B = 0. 

No quantitative theory has been developed for the crust motion in magnetars. The 
theory does not predict the geometry of the magnetosphere and its deformations, 
however, it is possible to predict what happens after a twist has been implanted. 
Part of the magnetosphere becomes filled with plasma and gradually dissipates the 
magnetic energy, creating a long-lived luminous corona. The theory of this activity 
can be tested with observations of the post-burst behavior of magnetars. Besides the 
nonthermal emission, bright hot spots are observed on magnetars after their bursting 
activity. The spots slowly shrink, on a timescale of months to years. Understanding 
this behavior may help reconstruct the geometry of the implanted twists and perhaps 
the crust motion that created them. 

The theory is developed in two steps. First, an electrodynamic model is formu- 
lated to describe the deformed magnetosphere and its evolution (§ 2). Second, basic 
properties of plasma filling the magnetosphere are evaluated (§§ 4 and 5). Observa- 
tional tests of the theory are discussed in §§ 3 and 6. 

In this contribution, we focus on the X-ray emission mechanism in magnetars. 
The related topics of their spindown behavior and optic al/IR/radio emission are dis- 
cussed elsewhere in this volume. 
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Twisted force-free magnetospheres were extensively studied in the context of the 
solar corona. In axial symmetry, they are described by the Grad-Shafranov equa- 
tion, and various force-free equilibria can be constructed numerically by solving this 
equation. A simple example is a self-similarly twisted dipole (Wolfson 1995). Pa- 
van et al. (2009) studied multipolar self-similar solutions. Realistic non-self-similar 
equilibria do not admit a simple description and are less explored, except for the 
case of weak twists (Wolfson & Low 1992; Beloborodov 2009, hereafter B09). 

The twisted magnetosphere stores additional energy - the magnetic energy of the 
electric currents maintaining V x B. The magnetosphere tends to dissipate the stored 
energy and untwists with time. Apart from the rare flares, the evolution proceeds 
slowly through a sequence of force-free equilibria. A key agent in this evolution is 
the moderate electric field En parallel to B. It has three important functions (BT07): 
(1) maintains the current flow, (2) regulates the dissipation rate Enj and the observed 
luminosity, and (3) determines the evolution of B in the untwisting magnetosphere. 

This section describes the untwisting theory for axisymmetric magnetospheres. 
A simple electrodynamic equation governs the evolution of B. 



2.1 Evolution equation for axisymmetric twist 

Any axisymmetric magnetic field B may be expressed in the following form, 

Vf x e<j> 

27Trsin0 Y Y 

Here r, 9,<p are spherical coordinates, is the angle of rotation about the axis of 
symmetry, and is the unit vector in the (j) direction. For example, a dipole mag- 
netosphere with dipole moment pL has f(r, 9) = 2%\i sin 2 9/r and B^ = 0. Function 
f{r, 9) has the meaning of the magnetic flux through the circle < <p < 2n of con- 
stant r and 9. Note that / is constant along any field line (B • Vf = 0) and by 
symmetry / is also constant on any 2D "flux surface" formed by rotating a field line 
about the symmetry axis. Any axisymmetric magnetosphere can be thought of as a 
set of nested flux surfaces labeled by /. Note that / = on the axis. 

If Bq = then the axisymmetric magnetosphere is potential, V x B = 0, and 
carries no electric currents. Any axisymmetric starquake is, in essence, a latitude- 
dependent rotation of the crust. It inevitably imparts a twist into the magnetosphere. 
Its amplitude is measured by the azimuthal extension of field lines outside the star, 

vCfl = te-fr= l QB i^ ( 2 ) 

Jp B rsmf? 



where the intergal is taken along the field line between its footprints P and Q on the 
star surface. The appearance of B§ ^ implies a non-zero circulation of B along 
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Fig. 1 Field line C\ with 
footprints Pi , Q\ on a flux 
surface f\ and its neighbor 
field line Ct with footprints 
Pi, Ql on a A ux surface ft- 
Closed contour C is formed 
by connecting Pi with P2 and 
Ql with Q2 along the star 
surface. 




IC, 



circles of constant r, 9. The circulation must be sustained by an electric current 
through the circle, according to Stokes' theorem. 

Current can flow only along B and must be maintained by a longitudinal voltage, 

rQ 

& e = J E ]{ dl, (3) 

where E\\ is the electric field component parallel to B. Note that both footprints P, 
Q of the field line are "grounded" in a good conductor — the neutron star — and 
the electrostatic voltage between P and Q is negligible. Voltage <t> e is inductive, not 
electrostatic. It is directly related to the rate of untwisting of the magnetic field, 
/ Endl = c~ l [ (dA\\/dt)dl where A is the magnetic vector potential (BT07). The 
gradual dissipation of the magnetospheric twist is, in essence, Ohmic dissipation in 
the electric circuit. The dissipation rate is proportional to <P e . 

How does the twist y/(/) evolve as a result of Ohmic dissipation? The electro- 
dynamic equation d B/dt = — cV x E states that electric field determines the evolu- 
tion of magnetic field. Its integral version (Faraday's law) can be used to derive the 
evolution equation for y/. Consider two infinitesimally close field lines C\ and C2 
on different flux surfaces f\ and /? = fi + Sf. A closed contour C can be formed 
by connecting the footprints along the stellar surface: Pi with P2 and Q\ with Qo 
(Fig. 1). Let 8<P be the magnetic flux through contour C. The Faraday's law states 

= _ / E dl,= [ Endl- f Endl, (4) 
c dt Jc Jc 2 JCi 

where we neglected the infinitesimal contribution from curves P1P2 and Q1Q2 to the 
line integral § c . Only contributes to flux 8<P. This can be seen if we continuously 
deform C\ and C2 so that their points slide at constant r, 9 in the -direction until 
both C\ and C2 end up in one meridional plane = 0. In this deformation, C\ and 
C2 slide along their flux surfaces f\ and f2, and hence 8<P remains unchanged by 
the deformation (magnetic flux through any part of a flux surface is zero). Only B§ 
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contributes to the magnetic flux through the deformed contour in the meriodinal 
plane. One can show that 80 is related to the twist angle \jf(f) by 

8* = ^. (5) 

This geometrical identity gives a quick way to derive the evolution equation for 
\lf(f,t). Combining equations Q, (0J, and © we immediately obtairQ 



= 2nc^-. (6) 
of 



Here we used d/dt = (d/dt)f taking into account that contour C remains on the 
same flux surfaces f,f + 8f as the magnetosphere untwists (the flux surfaces can 
change their shapes in the evolving magnetosphere). 

Equation (|6]l describes the twist evolution assuming that the footprints of field 
lines are frozen in the static crust of the star. The untwisting of field lines is due 
to their "slipping" in the magnetosphere, which changes the connectivity between 
the footprints so that y/ is reduced. It is straightforward to generalize the evolution 
equation for magnetospheres with moving footprints, 

-^-=2Kc-^ + (o{f,t) 1 (7) 

where (0 = d(j>g/dt — dfypjdt is the differential angular velocity of the moving crust. 



2.2 j -bundle 

Equation (O has an interesting feature: the evolution of y/ on a given flux surface / 
is governed by the transverse gradient of the voltage d<P e /df rather than its value 
<£>,,(/). This leads to a paradoxical result: if 3> e (/) = const i.e. the voltage is the 
same for all field lines, then the twist must "freeze" — the magnetic configuration 
does not change with time. How can this be reconciled with the continuing Ohmic 
dissipation in the electric circuit <P e j ^ 0? The circuit is ought to dissipate the mag- 
netic (twist) energy, and the unchanged B appears to violate energy conservation. 

The paradox is resolved if one notices that <t> e cannot be the same for all field 
lines. In particular, field lines / > fy are confined inside the star where En ~ 0. 
To fix ideas, suppose the star is a perfect conductor, which implies <P e (/) = for 
/ > /fi. Only field lines extending beyond the stellar surface, i.e. forming the mag- 
netosphere, have <P e 0, which implies a gradient d<P e / df^Q near the flux surface 
/ = fa. Then the twist evolution equation © implies that dy//dt < is maintained 



2 This evolution equation was derived in B09 for weak twists y/ <C 1 . As the derivation here shows, 
the same equation is also valid for strong twists. Equation (17) in B09 is not accurate for strong 
twists, it neglects the motion of flux surfaces in the process of untwisting. 
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near f = fit. Effectively, the magnetospheric electric currents are "sucked" into the 
star through the flux surface fa. This behavior is demonstrated by the explicit solu- 
tions to the evolution equation (|6]l in B09. 

There is another nontrivial feature of the untwisting process. Under normal con- 
ditions, plasma density in the closed magnetospherqj is not sufficient to conduct 
interesting electric currents. Plasma in the twisted magnetosphere is supplied by a 
continual discharge, which has a threshold voltage <P e = "V{f) ~ 10 9 V (BT07). The 
voltage is well described by a step function (see § 3 for details), 



The threshold nature of the discharge has a rather peculiar implication. Two distinct 
regions immediately form in an untwisting magnetosphere: "cavity" near the flux 
surface /r and "j-bundle" where j ^ (see B09 and Fig. 2). The two regions are 
separated by a sharp boundary along a flux surface /*. The untwisting process is the 
slow motion of the boundary /* that gradually erases the remaining electric currents. 
The erased magnetospheric currents end up closed below the stellar surface. 

If the discharge voltage (/) = const, the twist freezes in the j-bundle / < /* and 
passively waits for the front /* to come. Realistically, ^ const and d'f jdf > 
is possible. Then equation (O gives dy//dt > 0, i.e. the twist grows inside the j- 
bundle while waiting for the front to come. In all cases, the magnetic energy of the 
shrinking j-bundle decreases with rate equal to the rate of Ohmic dissipation. 

The behavior of untwisting magnetospheres is illustrated in B09 by solving the 
electrodynamic equation © for a concrete configuration: a centered dipole with a 
weak twist y <1. The "cavity + j -bundle" structure forms immediately, and then the 
cavity gradually expands, while the twist angle y/ inside the j-bundle grows linearly 
with time. This growth leads to i/a > 1 where the model breaks. General consider- 
ations (e.g. Uzdensky 2002) show that the magnetosphere with y/ > y/ cr = &{X) 
must lose equilibrium and restructure itself. Two possible types of restructuring 
are known: (1) Partial opening of field lines, which can lead to a configuration 
with lower energy (Wolfson & Low 1992). In an axisymmetric magnetosphere, this 
change may occur without breaking the symmetry. (2) Kink instability, which breaks 
the axial symmetry. Then an asymmetric plasmoid must be ejected by the magneto- 
sphere and the twist amplitude must be reduced. The exact threshold for the kink, 
the energy of the ejected plasmoid, and the new equilibrium configuration follow- 
ing the kink are unknown. One can expect that the amplitude of the growing twist is 
regulated by the instability so that it remains close to the critical value yi ~ 1 . 

The j-bundle and its footprints can be a bright source of radiation, with luminosity 
equal to the rate of Ohmic dissipation in the bundle. If the footprint of the j-bundle 
experiences heating by the bombarding magnetospheric particles (BT07) then /* 
defines the boundary of a hot spot on the star 6 < 0*. The area of this spot is given 



3 The usual estimate for plasma density is pcs/s where pcj = —£2 ■ B/2jtc is the charge density 
induced by rotation of the star with an angular velocity Q. (Goldreich & Julian 1969). 




■y j>o 

o j = 



(8) 




~ 5 lg(cj/BR) -2 0^ 1 

Fig. 2 Snapshot of the magnetospheric evolution. A global twist with a moderate y = 0.2 was 
implanted into the magnetosphere at / = 0, and the snapshot shows the magnetosphere at / ~ 1 yr. 
Details of the calculations are described in B09. The plane of the figure is the poloidal cross section 
of the magnetosphere. The black curves are the poloidal magnetic field lines; they are symmetric 
about the vertical axis and the equatorial plane. Left panel: current density. The cavity is shown 
in white (j = 0); its boundary /* is marked by the magenta curve. The cavity expands with time, 
and the j-bundle shrinks toward the vertical axis. Right panel: twist amplitude y/ at the same time. 
The neutron star is shown by the black circle. 



by 

A w 7r(i?sin0*) 2 = nR 2 u^. (9) 

Here we used the dimensionless variable u = f/ fu and approximated the poloidal 
structure of the magnetosphere by the dipole field0 As the cavity expands in an 
untwisting magnetosphere (Fig. 2), the spot 6 < 0* shrinks. 

The simplest model with "¥(u) = const gives the following estimate for the j- 
bundle luminosity (B09), 

L = Vh « ^Lyfuln 1.3 X 10 M V g B l4 R 6 y/ (^)' erg s" 1 , (10) 

where % = Y V and B\4 = B po i e /10 14 G. The evolution timescale of the lumi- 
nosity is given by 

t ev — 



dL/dt cRf 9 6r V0.1/ y 



Similar calculations have been done for the more detailed model where voltage "f 
is different on different field lines (B09). The theoretical expectations are compared 
with observations in the next section. 



4 Dipole field has f /2k = /J sin 2 6/r = /i/R max , where fl is the magnetic dipole moment. Then 
u = R/S max = sin 2 0] where d\ is the polar angle of the field-line footprint on the stellar surface. 
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3 Transient magnetars 

Recent observations revealed that many magnetars spend a significant fraction of 
time in a quiescent state with a low luminosity L ~ 10 33 — 10 34 erg s _1 and remain 
unknown until they produce an X-ray outburst. The post-outburst decay of the lumi- 
nosity back to the quiescent level is monitored in these "transient magnetars." The 
data can be compared with the theoretical model of untwisting magnetosphere. 



3.1 Magnetospheric activity or deep crustal heating? 

The mechanism of the magnetar bursts (e.g. Thompson & Duncan 1995; Woods & 
Thompson 2006) is not settled yet. It may be associated with starquakes. It may also 
be associated with sudden reconnection events in the magnetosphere that is slowly 
deformed by the plastic crustal flow (Lyutikov 2003). In either case, a strong and 
brief magnetospheric dissipation is invoked to explain the burst emission. 

Following the burst, a bright afterglow is observed, which continues to decay 
for months to years. The blackbody area of this emission is much smaller than the 
surface area of the neutron star, and it can be associated with the hot footprint of the 
j-bundle in the untwisting magnetosphere. This model relies on the external heating 
of the stellar surface by the energetic particles accelerated in the j-bundle. 

An alternative to external heating was considered by Lyubarsky et al. (2002). 
Their model assumed that during the burst a significant heat is deposited inside a 
deep region of the stellar crust. The subsequent passive cooling of this region due to 
heat diffusion through the crust can produce a bright thermal afterglow. Its expected 
spectrum is quasi-blackbody. 

Observations of several transient magnetars in the past few years give enough 
information to judge which of the two heating mechanisms — external or internal 
— is responsible for the long-lived afterglow. Two features point to external heating. 

• Nonthermal X-ray spectrum. The X-ray emission shows significant deviations 
from blackbody (BB) and is usually fitted by BB + power-law model. A high- 
energy X-ray component was discovered in the transient magnetar IE 1547.0- 
5408 (Enoto et al. 2010a). Its luminosity is comparable to the luminosity of the 
blackbody component, and the overall spectrum of IE 1547.0-5408 was similar 
to the spectra of magnetars with persistent coronal activity. The energy depen- 
dence of pulse profiles resembled that of persistent emission in 4U 0142+61 and 
1RXS J 170849-4009 10 (den Hartog et al. 2008a,b), suggesting a similar magne- 
tospheric emission mechanism. In another transient magnetar SGR 0501+4516, 
bright 100-keV emission was reported at t = 4 d after the outburst (Enoto et al. 
2010b), with a similar spectrum. 

• Shrinking area of the hot spot. The emission area of the blackbody component 
A decreases following the outburst together with its luminosity L. The transient 
blackbody emission evolves along tracks in the A-L plane as seen in Figure 3. 
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Fig. 3 Area A and luminosity L 
of the hot spots observed on five 
transient magnetars at different 
times following their outbursts; 
see text for references and dis- 
cussion of the individual ob- 
jects. One can see from the fig- 
ure that the spot area decreases 
together with its luminosity. The 
j -bundle model is shown by the 
strip between the two lines. It 
predicts the relation L = 1.3 x 



l^BuRfTgvtfi erg s" 1 . 
BuR^Tgij/ = 1 on the lower 
line and 20 on the upper line. 



10 11 
A [cm 2 ] 



The tracks are inconsistent with the model of impulsively heated deep crust, as 
cooling due to heat diffusion would give a hot spot that tends to spread with time, 
not shrink. The shrinking spot is consistent with the untwisting magnetosphere 
model (§ 2.2). The theoretical relation between A and L (obtained from eqs. [9] 
and[l0] see caption to Fig. 3) appears to agree with the data. Caution should be 
taken though when interpreting the details in Figure 3, as there may be signifi- 
cant systematic errors in observed A — it depends on the details of spectral fits. 
Besides, the distance d is poorly known for many magnetars. Note, however, that 
the shape of the observed track in the A-L plane does not depend on d, because 
both A and L scale as d 2 . 

These observations do not rule out deep crustal heating during the burst. The internal 
heating may dominate the thermal afterglow during first few days following the 
burst. After a few days, the observed light curve significantly flattens (Woods & 
Thompson 2006), suggesting a possible transition from internal to external heating. 



3.2 Individual objects 

XTE J1810-197 is a canonical transient magnetar. It has spin period P = 5.54 s and 
magnetic dipole moment fi w 1.5 x 10 32 G cm 3 , which corresponds to the surface 
field fi«3x 10 14 G. An X -ray outburst was detected from this object in January 
2004 (Ibrahim et al. 2004) and then its luminosity approximately followed an ex- 
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Fig. 4 Comparison of the model 
with the observed evolution of 
the area A and luminosity L of 
the hot spot after the outburst 
in XTE J1810-197. The data 
(open squares) are from Got- 
thelf & Halpern (2007). Dashed 
line shows the object luminos- 
ity in quiescence. Solid curves 
show the proposed theoretical 
model. The starquake occurred 
in the region u < uq = 0.15 
with amplitude y/o = 0.5. The 
discharge voltage is assumed to 
drop linearly from 5.5 GeV at 
u = 0.15 to 1 GeV at u = 0. 
(From B09.) 



2xl0 7 4xl0 7 6xl0 7 8xl0 7 10 8 

t [s] 



ponential decay on a timescale of 233 days for 3 years (Gotthelf & Halpern 2007). 
All data points for this object in Figure 3 were taken during the afterglow of this 
outburst. Remarkably, the source switched on as a bright radio pulsar (Camilo et 
al. 2007). The radio observations provided an estimate for the distance d = 3.3 kpc 
(which is used in Fig. 3) and accurate measurements of the spindown rate. In a few 
years, the object returned to its quiescent (pre-outburst) state: the radio pulsations 
and the spindown enhancement gradually disappeared, and the X-ray luminosity 
decreased to its quiescent level. 

These observations show that the magnetosphere of XTE J18 10-197 changed 
in the outburst. Its footpoints must have moved, imparting a twist in the magneto- 
sphere. The data give significant hints about the twist geometry and evolution: 

(1) The change in spindown rate suggests that the open field-line bundle was 
strongly affected, and hence the twist reached large amplitudes \j/ <; 1 near the mag- 
netic dipole axis. Recall also that y/ is limited to i// cl - = 0(1) by the MHD instability. 

(2) The luminosity was observed to decrease from ~ 10 35 erg/s to ~ 10 33 erg/s. 
The theoretically expected luminosity from a global twist with y/ ~ 1 would be 
much higher: L ~ 10 37 "/9 erg/s (eq. [TOb . We conclude that the twisted region was 
small: the current-carrying field lines formed a narrow bundle emerging from a small 
spot on the star surface. 
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(3) The spot was discovered: a hot blackbody component with a small emission 
area was found in the X-ray spectrum following the outburst (Gotthelf & Halpern 
2007; Perna & Gotthelf 2008; Bernardini et al. 2009). Its emission area shrank with 
time until the spot became undetectable (Fig. 4). 

(4) The X-ray and radio pulse profiles had almost simultaneous peaks, consistent 
with the X-ray emitting spot being on the magnetic dipole axis (Camilo et al. 2007). 

This behavior is consistent with the untwisting magnetosphere model described 
in § 2. The theoretical estimates for A, L, and f ev (eqs. [9] [TOj [TTJ are in remarkable 
agreement with the data. Figure 4 compares the detailed model with observations of 
XTE J18 10-197. The data constrain the magnetospheric voltage 'f to be between 1 
and 6 GeV, which is close to its expected theoretical value (§ 4). 

IE 1547.0-5408 is known as the magnetar with the shortest spin period P = 2 s. 
Its estimated magnetic field is B ~ 2 x 10 14 G. The object shares many features of 
XTE J1810-197. A similar dynamic hot spot was detected (Halpern et al. 2008). The 
data points for IE 1547.0-5408 in Figure 3 assume d = 4 kpc following Tiengo et 
al. (2010). One data point shows a pre-outburst state and two points show the decay 
after the outburst in 2007. The spot size was decreasing, however later the lumi- 
nosity and the spot grew again. In October 2008 new burst activity was detected, 
the emission area quickly increased to A 10 12 cm 2 and remained approximately 
constant for 20 days (Israel et al. 2010). A similar luminosity and area were ob- 
served after the outburst in January 2009 (Enoto et al. 2010a); this gave another data 
point in Figure 3. The activity level of IE 1547.0-5408 is clearly higher than that 
of XTE J1810-197. Apparently, many deformation events occurred in its magneto- 
sphere in the past few years, and it may be hard to see the long-timescale untwisting 
behavior because it is interrupted by new bursting activity. 

CXOU J164710.2-455216 is a transient magnetar with spin period P = 10.6 s 
and magnetic field B~9x 10 13 G. Its outburst in September 2006 was accompanied 
by the increase of persistent luminosity from ~ 10 33 erg s~ x to ~ 10 35 erg s _1 
(Muno et al. 2007). The subsequent slow decay of emission was observed for about 
150 days (Woods et al. 2010). All eight data points for this object in Figure 3 are 
based on these observations; distance d = 5 kpc was assumed following Muno et al. 
(2007). The decay timescale is comparable to that in XTE J1810-197. 

SGR 0501+4516 was discovered due to its bursting activity in August 2008 and 
then its post-burst behavior was observed (Rea et al. 2009). The results of these 
observations are shown by the five data points in Figure 3, assuming d = 5 kpc. The 
object has spin period P = 5.7 s and magnetic field fi~2x 10 14 G. Its post-outburst 
track in Figure 3 is consistent with a shrinking hot spot on the star. The track is close 
to that of XTE J1810-197. The luminosity decayed on a timescale of ~ 24 days, 
about 10 times faster than in XTE J18 10-197. This behavior is consistent with the j- 
bundle model if the voltage ~f in SGR 0501+4516 is ~ 3 times higher and the twist 
amplitude y is ~ 3 times smaller than the corresponding parameters in XTE J 1 8 1 0- 
197. This difference in parameters may not affect much the track in the A-L plane 
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since L=< y¥ (eq.fTOli. At the same time, it explains the short evolution timescale 
since f ev oc y- 1 f (eq.fTTT). 

SGR 0418+5729 was discovered due to its outburst in 2009 and its post-outburst 
behavior was followed for 160 days (Esposito et al. 2010). Its track in the A-L plane 
is shown in Figure 3 assuming a distance d = 5 kpc. The object has spin period 
P = 9.1 s and magnetic field fi£3x 10 13 G. The low magnetic field (by magnetar 
standards) makes this object special, which could explain the offset of its track from 
the tracks of other transient magnetars. Note also that the track location would be 
different for different d. A large distance d ~ 15 kpc would imply ~ 10 times larger 
L and A, and would place the track of SGR 0418+5729 into the strip between the 
two lines in Figure 3. The timescale of the post-burst decay is a few months (the 
luminosity decreased by a factor of 6 in four months). 

SGR 1627-41 is another transient magnetar whose spin period and magnetic field 
are close to those of IE 1547.0-5408. The object displayed bursting activity in 1998 
and then its luminosity gradually decayed toward the quiescent state with no new 
bursts detected until 2008 (Kouveliotou et al. 2003; Mereghetti et al. 2006; Esposito 
et al. 2008). The luminosity decayed in <~ 3 — 5-years, with a light curve similar 
to the post-outburst behavior of XTE J1810-197. The published spectral analysis 
provides almost no information on the blackbody emission area in this object except 
for one XMM observation when the source luminosity was only 4 x 10 33 erg _1 
(Mereghetti et al. 2006). 

A large fraction of observed magnetars are persistent bright sources with lumi- 
nosities L ~ 10 35 erg s _1 . The 100-keV emission component was observed in sev- 
eral of these "persistent magnetars," indicating quasi-steady magnetospheric activity 
(Kuiperet al. 2008). Their blackbody components have emission areas A <; 10 12 cm 2 
(e.g. Perna et al. 2001), significantly larger than in transient magnetars with simi- 
lar luminosities L ~ 10 35 erg -1 . The reason for this difference is presently unclear. 
Quasi-steady magnetospheric emission is expected if the twist decay timescale f ev is 
longer than the timescale of the crust motion imparting the twist. Note that f ev may 
be long if the voltage along the twisted field lines is somewhat below 10 9 V and/or 
the magnetic flux in the j-bundle is large (i.e. is large, cf. eq.QT). Moreover, it is 
possible that the j-bundle is persistently deformed by the plastic flow of the crust, 
and a quasi-steady twist is maintained near the threshold of MHD stability, l// ~ 1 . 

The luminosity of persistent magnetars L ~ 10 35 erg s (e.g. Durant & van Kerk- 
wijk 2006) is consistent with a twist created by rotation of a small cap on the star, 
with area A Ss 10 12 cm (see eq. [Tol and B09). Apparently, the footprints of the j- 
bundle in persistent magnetars are not strongly heated, as no hot spots with small A 
were detected in these objects. 

Persistent magnetars also display sudden increases and gradual decays of lumi- 
nosity. For example, the outburst of SGR 1900+14 in 1998 led to a rise in its persis- 
tent luminosity, and in a few months it decayed back to the pre-outburst level (Woods 
et al. 2001). Similar behavior was observed in AXPs, for instance in IE 2259+586 
(Woods et al. 2004) and 4U 0142+61 (Gavriil et al. 2009; Gonzalez et al. 2010). The 
rise in luminosity and its subsequent decay was reported for IE 1048.15937 (Dib et 
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al. 2009). The transient behavior in these objects is not necessarily associated with 
the same j-bundle that creates their persistent activity. The additional energy may 
be released on different field lines and after their untwisting the object returns to 
the pre-outburst persistent state. The timescale for this additional activity depends 
on the voltage, twist amplitude, and magnetic flux of the transiently deformed field 
lines. The variations in f ev may be similar to those in transient magnetars - at least 
a factor of 10. 

Observations of PSR J 1846-0258 suggest that the transient activity of the closed 
magnetosphere may also occur in rotation-powered pulsars (Gavriil et al. 2008). 



4 Pair creation 

Magnetospheric activity manifests itself through relativistic particles that can bom- 
bard the stellar surface and also generate nonthermal emission. We now discuss the 
origin of these particles, their energies, number, and circulation in the magneto- 
sphere. 

Particles are accelerated along the magnetic field lines by the longitudinal voltage 
<P e = J Ejjdl. The particles can resonantly scatter the stellar photons of energy hco 
once they reach the Lorentz factor 



which requires <P e ~ 10 9 /?i5 (Tio/lO keV)~' V. The scattered photons are boosted 
in energy by the factor of ~ y 2 . These high-energy photons convert to e ± pairs off 
the magnetic field, creating more particles. 

A similar process of e ± creation operates in the polar-cap discharge of ordinary 
pulsars, but in a different mode. In ordinary pulsars, the high-energy photons convert 
to e ± with a delay. The scattered photon initially moves nearly parallel to B and 
converts to e only when it propagates a sufficient distance where its angle Qy with 
respect to B increases so that the threshold condition for conversion is satisfied. This 
delay leads to the large unscreened voltage in the polar-cap models for pulsars. 

In contrast, the magnetic field of magnetars is so strong that pair creation can 
occur immediately following resonant scattering (BT07). The energy of the scattered 
photon is related to its angle Qy by 



where Eg = (2B/Bq + \) l l 2 m e c 2 is the energy of the first Landau level and Bq — 
m 2 ^ /he « 4.4 x 10 13 G. The scattered photon may immediately be above the 
threshold for conversion, E > E^ = 2m e c 2 / sin 9y, if B > 4-Bq. This suggests that e ± 
discharge in magnetars screens En more efficiently and buffers the voltage growth. 




(12) 
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Pair creation on field lines with apexes R max ^ 2R 

The discharge on closed field lines can be explored using a direct numerical ex- 
periment where plasma is represented by a large number of individual particles. 
The existing numerical simulations (BT07) describe the discharge on field lines that 
extend to a moderate radius R max 2R, where R is the radius of the neutron star. 
The magnetic field is ultrastrong everywhere along such field lines, B 3> Bq, and 
resonant scattering events may be effectively treated as events of pair creation - a 
significant fraction of scattered photons immediately convert to e ± . 

The simulations demonstrate that voltage and pair creation self-organize so that 
a particle scatters on average ~ 1 photon as it travels through the electric circuit, 
maintaining the near-critical multiplicity of pair creation ~ 1. This criticality 
condition requires voltage V ~ 10 9 V, which accelerates e ± particles to Lorentz 
factors y ~ 10 3 . The electric circuit operates as a global discharge, in the sense 
that the accelerating voltage is distributed along the entire field line between its 
footprints on the star. It is quite different from the localized "gap" that is often 
pictured above polar caps in pulsars. 

The discharge fluctuates on the light-crossing timescale ~ R/c and persists in 
the state of self-organized criticality. The behavior of the circuit resembles a con- 
tinually repeating lightning: voltage between the footprints of the field line quasi- 
periodically builds up and discharges through enhanced production of charges. The 
average plasma density in the circuit n is close to the minimum density n m ; n = j/ec, 
as required by the criticality condition . // = n/n m [ n ~ 1. 

As discussed in § 2.2, the currents tend to be quickly erased on field lines with 
small /? max - It may be that this part of the magnetospheres of observed magnetars 
is potential, V x B = 0, with no discharge. Then the observed activity is associated 
with currents on field lines with large 7? max , i.e. extending far from the star. 

Pair creation on field lines with apexes /? max > R 

The discharge on extended field lines is likely to have a similar threshold voltage 
f ~ 10 9 V, because the conversion of upscattered photons to e ± is efficient near 
the footprints where B ^> Bq. In this zone, particles are able to resonantly scatter 
soft X-rays once they are accelerated to y ~ 10 3 (cf. eq.fT2l. which requires voltage 
~ 10 9 V. Further growth of voltage should be stopped as it would cause the excessive 
creation of e ± moving in both directions, toward and away from the star, leading to 
efficient screening of E\\ . 

The discharge is expected to occur near the star, and some of the created particles 
outflow to r R along the extended field lines. The rate of resonant scattering by a 
relativistic particle increases as it moves from B 3> Bq to B Sa Bq. The particle scat- 
ters many more photons, because the resonance condition shifts toward photons of 
lower energy ha^^ <x B whose number density is larger. Note also that the effective 
cross section for resonant scattering, <7 res = 2n 2 r e c/(Oi &s increases as B l . Photons 
scattered in the region B ^ 10 12 G convert to e ± , which scatter more photons. 

In essence, the particles outflowing from the discharge zone lose energy to scat- 
tering as they enter weaker fields, and this energy is transformed to new generations 
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of e . A similar cascade above polar caps of pulsars was studied recently by Medin 
& Lai (2010). As a result, the e ± multiplicity of the outflow increases from ^ ~ 1 
to 3> 1 It means no charge starvation in the outer corona — there are plenty of 
charges to conduct the current demanded by the twisted magnetic field. 



5 Plasma circulation in the magnetosphere 

Dynamics of the e ± flow in the outer corona (r R) is strongly influenced by res- 
onant scattering, which exerts a strong force & on the particles along the magnetic 
field lines. & vanishes only if the particle has the "saturation momentum" p+ such 
that the radiation flux measured in the rest frame of the particle is perpendicular 
to B. In a weakly twisted magnetosphere with approximately dipole magnetic field 
exposed to central radiation is given b)Q 

2cos0 

^ = ^r (14) 

(momentum in units of m e c). The radiative force always pushes the particle toward 
p = p+. This effect may be measured by the dimensionless "drag coefficient," 



pm e c 2 



(15) 



Momentum p+ is a strong attractor in the sense that deviations p — p+ generate @ 3> 
1 in the outer corona. 

Consider a steady e outflow in the outer region where no new pairs are pro- 
duced. The outflow satisfies two conditions: (a) it is nearly neutral, n+ « «-, and 
(b) it carries the electric current j demanded by V x B. Electric field E\\ immediately 
reacts to violations of these conditions and enforces them. From (a) and (b) one can 
derive the relation between the velocities j3± of e ± (details are given in Beloborodov 
2010), 

0- 2 

Here "+" corresponds to positrons, which carry current j + > and "— " corresponds 
to electrons, which carry j- < 0; we assume a positive net current j = j + + for 
definiteness. The multiplicity is defined by ./# = (j + + \j-\)/ j- An electric field 
must be generated in the outflow to sustain the condition ( [ToT l against the radiative 
drag that tends to equalize j3_ and j8 + at j3*. This electric field is generated by a 
small deviation from neutrality, 8n = n + — «_ n + . 



5 A steady relativistic outflow without pair creation would have = const along the field lines. 
This follows from conservation of magnetic flux, charge, and particle number, which give n/B = 
const, j/B = const, and nj j — const along a field line. 

This expression is valid in the region where \-B r jB > {R/r) 2 . In this region, stellar radiation 
may be approximated as a central flow of photons, neglecting the angular size of the star ~ R/r. 
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Fig. 5 Lorentz factors y + (left panel) and y_ (right panel) in the two-fluid model of outer corona. 
In this example, the electric current is carried by the e outflow of a fixed multiplicity = 50. 
The plasma is injected at radius r = 2R and outflows along the magnetic field lines (white curves). 
The flow is illuminated by the star with temperature kT = 0.5 keV (magenta circle at the origin), 
which exerts radiative forces ^{j±) on the e ± plasma. The Lorentz factors y + and y_ change as 
the flow enters the drag-dominated region | & \ ~S> 1 . The region \S\ > 3 is shown by the thick black 
curve and shadowed in black. & < for positrons (y + > y+) an d & > for electrons (y_ < y*). 
The radiative drag stops the plasma in the equatorial plane outside ~ SR. A nearly dipole magnetic 
field (weakly twisted) with B po ie = 10 15 G is assumed in this example. R is the neutron-star radius. 



The two-fluid dynamics of the outflow is governed by equations, 

mec2 ^df = ^(y+) +eE b mec2< ^r = &(y-)- eE b ( 17 > 

where / is the length measured along the magnetic field line. Since y_ and y + are 
not independent — they are related by condition ([Tol l — it is sufficient to solve one 
dynamic equation, e.g. for y + (and use the dynamic equation for y_ to exclude En). 
Straightforward algebra gives, 

n^**± = *^ + fb-\ d -^=(^-) 2 ( Y -)\ (18) 



dl l+dj-/dy+ dy+ \.M +\ ) \r 

In the region of strong drag, \@\ 3> 1, one finds i^(y+) « — &(Y-) ~ — eEn. Here 
both e + and e _ components of the outflow have velocities that are "locked" by the 
balance of two strong forces — radiative and electric. 

To illustrate the drag effect on the e flow around magnetars, consider the fol- 
lowing simplified model. Suppose that the e ± plasma is injected near the star with 
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a given multiplicity, e.g. „# = 50, and a given high Lorentz factor, e.g. y + = 100. 
The corresponding y_ is determined by equation ( fT6l ). Suppose that the plasma is 
illuminated by the blackbody radiation of the star of temperature kT = 0.5 keV and 
neglect the radiation from the magnetosphere itself (see below). The steady-state 
solution for y± established in the outer corona is shown in Figure 5. The solution 
is not sensitive to the precise radius of e ± injection as long as it is small enough, 
before the plasma enters the drag-dominated region. The most significant electric 
field develops in the region where \&\ ^> 1. It is given by eE» « —JP(y + ), and the 
corresponding longitudinal voltage established in the outer corona is found by in- 
tegrating J^(y+) along the field line, e<P e ps — / ,^(y + )dl. Its typical value for the 
model in Figure 5 is ~ 10 7 V. Flows with lower j% develop stronger electric fields, 
however in all cases of interest (j#t ^> 1) the drag-induced voltage is below 10 9 V. 

The calculations shown in Figure 5 assume that the plasma is optically thin to 
resonant scattering. This is not so for real magnetars. TLK02 showed that the optical 
depth of a strongly twisted magnetosphere is comparable to unity. When the large 
pair multiplicity is taken into account, the estimate for the optical depth becomes 



The large optical depth has two implications. First, there is a lot of scattered radi- 
ation in the magnetosphere, which is quasi-isotropic. This increases the drag and 
reduces p+. Second, there is a self-shielding effect: the drag force & experienced 
by an electron (or positron) is reduced by the factor of T . These two effects will 
change y±(r, 9) shown in Figure 5. However, the main feature will persist: the out- 
flow remains to be drag-dominated in the equatorial region at r ~ (10 — 20)R. The 
initial momentum of the e flow ejected from the discharge zone to the outer corona 
is taken away by the radiative drag and the e pairs accumulate near the apexes of 
the closed magnetic field lines. The pairs annihilate there. 

The annihilation rate is given by N ann ~ 2„#(//e). Here / is the electric current 
through the annihilation region, / ~ (ycf! /4R^), \j/ is the twist amplitude, Ri is the 
inner boundary of the annihilation region, and jj. is the magnetic moment of the star. 
The corresponding annihilation luminosity is given by 



The proposed picture of the magnetar corona is schematically shown in Figure 6. 
The exact position R\ of the boundary between the two regions in Figure 6 has not 
been calculated and may be inside SR. Consider, for instance, the plasma flow along 
the field lines with R max — 5R. It experiences scattering with rate that is high enough 
to generate a large e ± multiplicity, however, not enough to bring the plasma to a stop 
at the equatorial plane. The opposite flows in the northern and southern hemispheres 
meet in the equatorial plane and try to penetrate each other. The two-stream instabil- 
ity then develops and a strong Langmuir turbulence is generated, which will inhibit 
the penetration. Effectively, the opposite flows collide in the equatorial plane and 



>1. 



(19) 




(20) 
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Fig. 6 Schematic picture of plasma circulation in the magnetosphere with surface B ~ 10 15 G. 
Two regions are indicated. (1) "Inner corona." Here e have a moderate multiplicity ~ 1. 
The particles do not stop in the equatorial plane. The electric field Ey ensures that electrons and 
positrons circulate in the opposite directions along the magnetic field lines, maintaining the electric 
current demanded by V x B. The particles are lost as they reach the footprints of the field line and 
continually replenished by pair creation. (2) "Outer corona" — extended field lines with S max 2> R. 
Electrons and positrons are created by the discharge near the star and some of them flow outward 
to the region of weaker B. Here resonant scattering enhances the pair multiplicity, ^ S> 1, and 
decelerates the outflow. The e particles stop at the apexes of magnetic field lines (blue region 
in the equatorial plane), accumulate, and annihilate there. The number fluxes of electrons and 
positrons toward the annihilation region differ by a small fraction ~ ^tft~ 1 , so that the outflow 
carries the required electric current j = (c/4jt)V x B. 

should stop there. This behavior contrasts with the inner corona where the counter- 
streaming e~ and e + never stop; their circulation is maintained by the induced elec- 
tric field. The outer corona is different because of its high pair multiplicity then 
the electric current is organized while both e~ and e + outflow from the star. 

The above analysis assumed a steady e outflow. Then the current is organized by 
locking the electrons and positrons at two different velocities /3_ and j3+. In reality, 
the two-fluid flow is prone to two-stream instability that generates strong Langmuir 
waves and broadens the momentum distribution of e . The robust generation of 
plasma waves gives a mechanism for low-frequency emission from magnetars (dif- 
ferent from drag-induced waves in normal pulsars, cf. Lyubarsky & Petrova 2000). 

The instability has been simulated numerically, and the results are described in 
detail elsewhere (Beloborodov 2010). The fluctuating electric field £| is so strong 
that it could stop a particle with Lorentz factor y + on a short scale, much shorter than 
the free path to resonant scattering. A strong anomalous resistivity could be expected 
in this situation. Surprisingly, the simulations indicate that the effective resistivity 
(and the corresponding voltage) remain low. A complicated time-dependent pattern 
is organized in the phase-space, which allows the charges to find small-resistance 
paths through the waves of E\\ and conduct the current at a low net voltage and a 
low dissipation rate. 
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6 Magnetospheric emission 

Various ideas have been proposed for the mechanism of hard X-ray emission from 
magnetars (Thompson & Beloborodov 2005; Heyl & Hernquist 2005; BT07; Baring 
& Harding 2007; Lyubarsky & Eichler 2008). 

One possible source is the transition layer between the corona and the star. The 
layer can reach temperatures kT <; 100 keV and produce bremsstrahlung emission. 

The other possible source is the extended corona itself. It generates nonthermal 
X-rays via resonant scattering of the thermal X-rays emitted by the neutron star. A 
simplest model assumes that the high-energy spectrum forms via single scattering 
by an optically thin flow that interacts with the stellar radiation. This model was de- 
veloped for ordinary pulsars (e.g. Kardashev, Mitrofanov, & Novikov 1984; Daugh- 
erty & Harding 1989; Sturner 1995; Lyubarsky & Petrova 2000) and also applied 
to magnetars (Baring & Harding 2007). It does not, however, explain the magnetar 
spectrum. Several works emphasized that radiative transfer with multiple resonant 
scattering in the closed magnetosphere has a strong impact on the observed X-ray 
spectrum (TLK02; Lyutikov & Gavriil 2006; Fernandez & Thompson 2007; Nobili, 
Turolla & Zane 2008; Rea et al. 2008; Pavan et al. 2009; Zane et al. 2009). The key 
unsettled problem is what parameters of the scattering plasma should be assumed in 
the transfer calculations. Resonant scattering depends on the Lorentz factors of the 
coronal particles and the direction of their motion. 

One popular assumption is that the corona is filled with counter-streaming pos- 
itive and negative charges with mildly relativistic velocities. The scattering in this 
picture can reproduce the 1-10 keV part of the magnetar spectrum (e.g. Rea et al. 
2008). However its theoretical basis is problematic. The keV photons are resonantly 
scattered at radii r ~ 10/? where B ~ 10 11 — 10 12 G. The plasma strongly interacts 
with the stellar radiation in this region, and the counter- streaming model needs £| 
to push charges of one sign toward the star against the radiative drag. This electric 
field accelerates the charges of the opposite sign away from the star. In the presence 
of e ± plasma (which is inevitable), no self-consistent solution exists for the mildly 
relativistic counter-streaming model. 

A different picture of plasma circulation is shown in Figure 6. We argued that in 
the outer corona both + and — charges must outflow from the star and annihilate in 
the equatorial plane of the magnetic dipole. A moderate electric field E\\ is generated 
to maintain a small difference Aj5/f5 ~ <C 1 between the average velocities 
of the positive and negative charges, so that the outflow carries the required electric 
current. This picture offers a self-consistent solution for the drag-dominated electric 
circuit, which needs to be tested against observations. The theoretical X-ray spec- 
trum can be found by solving for the transfer of the stellar radiation through the e ± 
outflow. Then one can compare the theoretical prediction with the data. 

The technical difficulty is that the outflow dynamics is coupled to the radiation 
field, and hence the radiative transfer must be solved together with the plasma mo- 
tion. This nonlinear problem is well defined and can be solved exactly (numerically), 
however this is a rather formidable task. As a first step, let us consider a simplified 
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Fig. 7 Magnetar spectrum 
formed by resonant scattering 
in the e outflow. Solid his- 
tograms show models A and B 
(see text). The line of sight is 
chosen at angle = 60° with 
respect to the magnetic axis. 
The star is assumed to emit 
the blackbody radiation with 
temperature kT = 0.5 keV. The 
result of single (first) scattering 
in Model B is shown by the 
dotted histogram. 



1Q 34 _ 



0.1 



10 



100 
E [keV] 



1000 



10* 



model which provides the Lorentz factor y(r,9) of the scattering medium, moti- 
vated by the results of § 5. Given y(r, 9) it is straightforward to calculate the transfer 
of the stellar radiation through the corona. The result will not depend on the opti- 
cal depth T as long as T > 1 (this is a special feature of resonant scattering). We 
use the Monte-Carlo method for the transfer calculations; it is described elsewhere 
(Beloborodov 2010). 

First, consider the outflow model in Figure 5 and suppose y{r,9) — y+(r, 0) 
(model A). The emerging spectrum after multiple scattering in this outflow is shown 
in Figure 7 for an inclination angle = n/3. The model is inconsistent, as Figure 5 
assumed that the drag force is created by the central thermal radiation field, ne- 
glecting the scattered radiation, which would be appropriate only if T <C 1. The 
central radiative drag leads to the artificially high y + near the axis where the central 
radiation is beamed along the outflow direction. 

In reality, T> 1 (even for moderate twists, see eq. [19), the scattered radiation 
fills the magnetosphere and contributes to the drag force, decelerating the outflow 
near the axis. This will lead to small y(r, 9) at all 9 in the outer corona and change 
the scattered radiation spectrum. To see the trend of this change, let us modify the 
drag force & in Figure 5. Let a be the angle between the photons and the outflow 
velocity. In the region where cos a > 0.5 for all photons emitted by the star, let us 
set cos a = 0.5 when calculating jF. This intrusion makes the radiation field less 
anisotropic so that it is never strongly beamed along the outflow. Re-calculating 
Y±(r, 9), we find that they are reduced compared with those in Figure 5. Now con- 
sider the scattering medium with the new Lorentz factor y = y + and calculate the 
radiative transfer through this medium (model B). The result is shown in Figure 7. 
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Toy models A and B allow one to make some preliminary conclusions. In both 
models, the spectrum has two distinct peaks, at a few keV and near 1 MeV, and a 
minimum near 10 keV. This shape appears to be a robust result of radiative transfer 
through the e ± outflow. The result is consistent with the observed spectra of mag- 
netars (e.g. Kuiper et al. 2008). The photon index F above 10 keV is close to the 
typical observed F ~ 1 . 

Multiple scattering is essential for the formation of the high-energy peak. For il- 
lustration, the dotted histogram in Figure 7 shows the spectrum obtained in model B 
after single scattering; its high-energy component is weak. The first scattering of 
the thermal keV photons occurs mainly at radii r ~ \QR where the outflow is slow. 
It generates quasi-isotropic radiation with a spectrum that is slightly modified from 
blackbody. A fraction of the scattered photons propagate back to the star and ap- 
proach head-on the relativistic outflow in the inner region. The second scattering 
for such photons dramatically boosts their energy, contributing to the high-energy 
component in the observed spectrum. This mechanism will also operate in the exact 
model where the accurate distribution of the Lorentz factor in the outflow is found 
consistently with the radiative transfer. In this model, the high-energy luminosity 
will be explicitly fed by the initial power of the outflow injected near the star (which 
is determined by the electric current and the discharge voltage). The radiative drag 
transforms this power to radiation. 

Similar spectra are generated when the outflow is confined to the j -bundle of 
moderate magnetic flux m* <; 0.1. The spectrum evolves as the j -bundle shrinks to 
m* -C 0.1. Note also that the high-energy component is anisotropic — it is prefer- 
entially beamed along the magnetic dipole axis. This may help explain the energy 
dependence of the pulse profiles in observed magnetars (den Hartog et al. 2008a,b). 

The pulse profiles of magnetars are expected to change as the magnetosphere 
untwists and the j-bundle shrinks. Previous discussion of pulse profiles assumed 
that the entire magnetosphere is twisted in a self-similar way and has a moderate 
optical depth t proportional to the amplitude of the global twist y (e.g. Fernandez 
& Thompson 2007). The calculations presented here suggest a different picture. The 
twists of magnetars quickly evolve to narrow j-bundles (and twists may always be 
confined to a small part of the magnetosphere). The plasma in the j-bundle is opaque 
to resonant scattering even for weak twists y/ ~ 0.1. The result of radiative transfer 
(and the corresponding pulse profile) does not directly depend on T as long as T 3> 1 . 
It may depend on X indirectly, because the deceleration of the e ± outflow depends 
on its optical depth. Further analysis is needed to explore the implications for the 
pulse profile and to test the model against observations. 

A robust result of the radiative transfer (especially in the more reasonable 
model B) is the sharp break in the spectrum at ps 2m e c 2 w 1 MeV. It is caused 
by photon absorption in the strong magnetic field. The emission of multi-MeV pho- 
tons occurs in the region where B ~ 10 13 G, and they see a large optical depth to 
absorption. Absorption is reduced for photons emitted along the magnetic axis, and 
the emerging high-energy spectrum depends on the line of sight. 

Note that somewhat smaller £br ~ 0.5 MeV could be possible for the emission 
produced by the transition layer between the corona and star (Thompson & Be- 
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loborodov 2005; BT07). In that model, may be controlled by the temperature of 
the layer rather than absorption in the magnetosphere. Observations of E\, r may help 
distinguish between the two mechanisms of high-energy emission. 



7 Conclusions 

Nonthermal X-ray emission and shrinking hot spots on magnetars are associated 
with magnetospheric activity — the release of the energy stored in the twisted closed 
magnetosphere. The rate of energy release is proportional to the voltage established 
along the magnetic field lines, which is regulated by the continual discharge to 
<P e ~ 10 9 V. This voltage has a pure inductive origin, directly related to the rate 
of untwisting of the magnetic field. For axisymmetric configurations, the untwisting 
evolution is described by the electrodynamic equation (0. 

The untwisting theory predicts the formation of shrinking hot spots on "transient 
magnetars" whose magnetospheres are temporarily activated and gradually relax 
back to the quiescent state. Shrinking hot spots were indeed reported in these objects 
(Fig. 5), and their evolution appears to agree with the theoretical expectations. 

Self-similar twists (Wolfson 1995; TLK02) do not form in magnetars, because 
the electric currents tend to be quickly removed from field lines with moderate apex 
radius R max . Currents have longest lifetime on field lines with R max ^ R, forming 
the extended j-bundle. An e ± outflow of a high multiplicity streams in the j-bundle 
and creates the outer corona around the neutron star. The outflow stops at the apexes 
of the closed magnetic field lines and the e ± pairs annihilate there. 

The e outflow is opaque to resonant scattering and impacts the observed X-ray 
spectrum. The simulations of radiative transfer through the outflow suggest a natural 
explanation for the hard X-ray component in magnetar spectra. 

The strong radiative drag and the imposed electric current lock the j-bundle 
plasma in a peculiar two-fluid state. The two-stream instability is then inevitable; it 
creates strong Langmuir oscillations, which can convert to escaping low-frequency 
radiation. This mechanism is discussed in more detail elsewhere (Beloborodov 
2010). 
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